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Abstract 

We study formation of molecules in primordial prestellar clumps and evaluate the 
line luminosities to assess detectability by next-generation facilities. If the initial 
H2 fraction is sufficiently high, HD becomes an important coolant in the clumps. 
The luminosity from such HD cooling clumps is lower than that from H2 cooling 
ones because of lower temperature (< lOOK). As for Li reactions, we include the 
three-body LiH formation approximately. The Li molecular fraction remains very 
low (< 10~^) throughout the evolution owing to the high dissociative reaction rate of 
LiH + H ^ Li + H2. LiH does not become an important coolant in any density range. 
The luminous emission lines from the prestellar cores include H2 rovibrational lines: 
1-0 Q(l), 1-0 0(3), 1-0 0(5), and pure rotational lines: 0-0 S(3), 0-0 S(4), 0-0 S(5). 
The next-generation facilities SPICA and JWST are able to detect H2 emission in 
a large pre-galactic cloud that forms metal-free stars at a high rate of ~ lO'^Moyr"^ 
at redshift z < 10. We also derive an analytical expression for the luminosity that 
reproduces the numerical results. 

Key words: cosmology: early universe — galaxies: high-redshift — infrared: 
galaxies — molecular processes — stars: formation 



1. Introduction 

Population HI (Pop HI) stars, which are formed from primordial pristine gas, have been 
postulated as the origin of metals of metal-deficient stars in the Galactic halo (Schwarzschild, 
Spitzer 1953). Recent detection of very metal-poor stars renewed an interest in Pop HI stars 
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(e.g. Christlieb et al. 2002). Pop III stars may also account for the presence of metals 
3.5 X 10"^Zo) in the high-z {z ~ 5) La absorption systems observed in quasar spectra (e.g., 
Songaila 2001). Additionally, as a source of ultraviolet photons that caused the reionization of 
the intergalactic medium at z ~ 17±5 suggested by WMAP data (Bennett et al. 2003), Pop III 
stars are attracting the attention of many authors (e.g., Ciardi et al. 2003). These observations 
clearly suggest the importance of roles played by Pop III stars in the early evolution of galaxies. 

Partly stimulated by the observations, theoretical studies on the formation of primordial 
objects in the early universe have made remarkable progress in recent years. Simultaneously, 
several next-generation large observational facilities are planned, one of whose goals is to detect 
star formation in the first galaxies. In this context, a variety of observational verifications for 
scenarios of the first object formation in the early universe has been proposed. 

Detectability of primordial pre-galactic clouds has been discussed by some authors. For 
example, Haiman, Spaans, and Quataert (2000) and Fardal et al.(2001) estimated the La lumi- 
nosity from collapsing pre-galactic clouds and discussed their detectability by the James Webb 
Space Telescope (JWST) ^. Omukai and Kitayama (2003) demonstrated that H2 lines, in addi- 
tion to La emission, were emitted in collapsing pre-galactic clouds and discussed detectability 
of H2 rotational lines by the Single Aperture Far-Infrared observatory (SAFIR) ^. H2 line emis- 
sion from the supershells produced in IGM by supernova (SN) blastwaves of Pop III stars can 
be detected by JWST (Ciardi, Ferrara 2001). Observational features of clusters of Pop III 
stars that have already formed have been also studied by many authors (e.g., Bromm et al. 
2001; Schaerer 2002; Venkatesan et al. 2003). Near-infrared background anisotropies produced 
by these clusters may have been detected by the Two Micron All Sky Survey (2MASS) (e.g., 
Salvaterra, Ferrara 2003; Kashlinsky et al. 2004). 

H2-line luminosity from individual Pop III star-forming regions has also been evaluated 
(Kamaya, Silk 2002; Ripamonti et al. 2002; Mizusawa et al. 2004, hereafter paper I). At first, 
the H2-line emission from the primordial cloud including a number of first star-forming regions 
was expected to be observable by next generation facilities (Kamaya, Silk 2002; Ripamonti et al. 
2002). Later, a flaw was found in the luminosity distance used in their estimations. With the 
correct luminosity distance, detectability of the H2 emission at the typical epoch for first star 
formation [z ~ 20) is extremely low even with the Space Infrared Telescope for Cosmology and 
Astrophysics (SPICA) ^, which is one of the most sensitive facilities planned in this wavelength 
band (paper I). 

Nonetheless, even in lower redshift (say, 2 ~ 5) primordial gas can be present in objects 
that are formed with delay (Scannapieco et al. 2003). Detection of primordial molecular 
emission from the first episode of star formation in such objects should be easier than in genuine 

^ http://www.jwst.nasa.gov/ 

^ http://safir.jpl.nasa.gov/whatls/index.asp 

^ http:/ /www. ir.isas.ac.jp/SPICA/ 
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first objects in the universe. In the case of star formation in the first objects, whose mass is 
small (~ IO^Mq), dense star- forming cores are formed with httle fragmentation (Bromm et al. 
1999,2002; Abel et al. 2000,2002). On the other hand, in large pre-galactic clouds, post-shock 
cooled sheets, which are expected to form after virialization, fragment into filamentary clouds 
(Yamada, Nishi 1998; Yamada 2002). After further contraction and fragmentation of these 
filaments, star-forming cores are produced. In such clouds, H2 and HD form abundantly owing 
to ionization accompanying the strong shock (Uehara, Inutsuka 2000; Nakamura, Umemura 
2002) and strong molecular emission is expected. Another trace molecule in primordial gas, 
LiH, can be an important coolant if all the Li become LiH in spite of small Li abundance (Lepp, 
ShuU 1984). Luminosities of these trace primordial molecular (i.e., HD and LiH) emission were 
calculated, and their detectability by ALMA was discussed by Kamaya and Silk (2003). In 
their estimation, however, they again used the incorrect luminosity distances as in Kamaya and 
Silk (2002). Moreover, their assumption that all the deuterium and lithium are in molecular 
form, in particular for Li, is not justified, as described later. 

In this paper, we investigate the formation of HD and LiH as well as H2 not also in 
collapsing star-forming cores but also in filamentary clouds, and estimate their line luminosities. 
Here, to assess the importance of HD and LiH lines as coolants, we include the deuterium and 
lithium chemical reactions in addition to the processes considered in paper I and calculate 
these molecular fractions and luminosities. This treatment improves that of Kamaya and 
Silk (2003), who adopted a rather crude assumption for HD and LiH fraction without solving 
chemical reactions. In particular, we include the Li three-body reaction in an approximate 
way, by which the complete conversion of Li into LiH has been envisaged (Stancil et al. 1996). 
Additionally, we consider the formation reaction of LiH 

Li + H2^LiH + H, (1) 

which is the inverse reaction of the dominant LiH dissociation reaction. Since their rate coeffi- 
cients are not well contrained, these reactions have not been included in studies on the first star 
formation. Considering possible variations of these rates, we study the effect of different rate 
coefficients. Furthermore, we also assess the effects of different initial conditions, which corre- 
spond to different sizes of parent clouds, while in paper I we only calculated a single case. As for 
the geometry of clouds, we consider not only spherical cores, which are assumed in paper I, but 
also filamentary clouds, which are suppposed to form in large pre-galactic clouds (^ IO^Mq). 
As in paper I, we evaluate the luminosities of primordial molecules both in the accretion and 
collapse phases, while related works by Kamaya and Silk (2002,2003) and Ripamonti et al. 
(2002) studied only the collapse phase. 

Using this model, we discuss the condition, i.e., formation redshift, star formation rate 
(SFR), and size of parent pre-galactic clouds of Pop HI stars that permit SPICA and JWST 
to detect primordial molecular emission. We found that those facilities can observe H2 line 
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emission in a pre-galactic cloud at 2; < 10 with SFR ~ 10^'"^ Mq/jt. LiH fraction is found to 
remain low (< 10"'^ of total Li nuclei) even in such high densities where the three-body reaction 
is important. This is because of the high rate of dissociation reaction 

LiH + H^Li + Ha. (2) 

The structure of this paper is as follows: In section 2, our calculation method and 
results for spherical star-forming cores are described. In section 3, modification in the method 
and results for filamentary clouds are presented. In section 4, we give discussions about (i) an 
analytical estimate of molecular line luminosities, which agrees well with numerical calculations 
in sections 2 and 3, and (ii) detectability of primordial star-forming regions by next-generation 
facilities, SPICA and JWST. Finally, we summarize this paper briefly in section 5. 

2. Evolution of a Spherical Star-Forming Core 

2.1. The Model 

In this paper, we use the same method as in paper I in calculating the dynamical and 
luminosity evolution of star-forming cores. Modification from the previous model is addition of 
deuterium and lithium chemical reactions and the cooling due to HD and LiH. 

2.1.1. Dynamics 

In our analysis, dynamics in two phases of star formation, i.e., the runaway collapse 
phase and the main accretion phase, corresponding to before and after the formation of a central 
protostar, are treated in different ways because of their vastly different dynamical natures. 

In the run-away collapse phase, the dynamics is described mathematically by a Larson- 
Penston type similarity solution (Larson 1969; Penston 1969). In this solution, the cloud can be 
divided into two regions: the central densest region and the lower density envelope. Since the 
central densest region collapses at roughly the free-fall time, = ^3n/32Gpc, which inversely 
depends on the density, the envelope with lower density remains almost unevolved. Thus this 
evolution is called the run-away collapse. The central region with a flat density distribution 
has a length scale about the local Jeans length, Aj = ^Jcln/Gpc, where Cg and pc are the sound 
speed and density at the center. On the other hand, the envelope is characterized by density 
decreasing with radius by the power-law p oc r~'^^^'^~"'\ where 7 = d\nP/d\np is the effective 
adiabatic coefficient (Larson 1969). 

Here, we model the central region as a sphere with radius -Rcore = «Aj/2, and collapsing 
at the free-fall rate with pressure modification: 

dpc^Jh_ /ON 

dt pte' ^ ' 

where a and (3 are correction factors of order unity owing to finite pressure gradient force. 
These values will be given below shortly. Considering the run-away nature of the collapse, we 
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assume that physical quantities evolve only when the mass shell is inside the central region, 
i.e. r < Rcore'- after the matter has left behind the central evolution and become part of the 
envelope, those quantities are kept constant until the end of this phase. With this assumption, 
we can construct the whole density and velocity structures by calculating only the evolution of 
central quantities. 

The evolution of primordial star-forming clouds in the run-away collapse phase is studied 
by way of one-dimensional (spherical symmetric) radiative hydrodynamics by Omukai and 
Nishi (1998). They found that actual dynamics is approximated by the Larson-Penston type 
similarity solution (Larson 1969; Penston 1969) with 7= 1.09. For the self-similar solution with 
7 = 1.09, a = 1.13 and (3 = 2.13, which we use hereafter. 

At the end of the collapse phase, a small protostar forms at the center. Although the 
protostar is initially small, it grows in mass by accreting the surrounding matter. This is the 
onset of the main accretion phase. In this phase, since the pressure gradient force is negligible 
compared with the gravitational force, we assume that the matter in the envelope infalls at the 
free-fall rate. The evolution of each mass shell is given by 
dv _ GMr 

It ~ ^ ' 

where v is the velocity of a free-falling mass shell, and is the mass contained inside the 
mass shell. The density of a free- falling mass shell, p, is described by p oc r~^/^. Using these 
equations, we calculate the evolution of each mass shell until the begining of H2 dissociation. 

2.1.2. Thermal Evolution 

We follow the evolution of the internal energy per unit mass 

' (6) 



7ad - 1 pmH 
by solving the energy equation 

dt dt\p) p ' 

where P is the pressure, /i is the mean molecular weight, and mn is the mass of the hydrogen 
nucleus. The total cooling rate per unit volume, Atot, is the sum of the following contributions: 

Atot = Ah2 + AhD + ALiH + Achcm, (8) 

where Ahj is the cooling rate by H2 line emission, Ahd by HD line emission and Alih by LiH line 
emission, and Achem is the net cooling rate by chemical reactions. The dominant contribution 
to Achcm is the cooling by H2 dissociation. Although we have also included the cooling by 
continuum emission as in paper I, it is always negligible in comparison with other terms. The 
H2 cooling is treated in the same way as in paper I, where we considered effects of radiation 
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Tabic 1. Additional Li Chemical Reactions. 



Numbers 


Reactions 




Rate CoefBcients 


1 


Li - 


h H2 ^ LiH 


+ H 


fci = 2.0 X 10-"(QLiH/QH2)exp(-2.28 x IOVT) cm^ f 


2 


Li - 


h 2H LiH 


+ H 


fcs = 2.5 X 10-'^^ /T cm^ 


3 


Li - 


h H + H2 ^ 


LiH H 


- Hs fcg = 4.1 X cm^ 



Table 1. In addition to tiic Li reactions except photorcactions in Table 1 of Stancil, Lepp, and Dalgarno (1996), we included 
the three reactions above. 

t Qh2 ^^'^ QhiH a^ro partition functions of H2 and LiH from Sauval and Tatum (1984). 



trapping for optically thick lines by the escape probability method. In addition to the above 
processes considered in paper I, here we include the HD and LiH cooling. For the collisional 
de-excitation rate of HD, we use the rate coefficients of Flower et al. (2000). We take account 
of the first eight rotational levels of HD. In calculating the emissivity of HD lines, we use the 
same method as for H2 lines. The cooling rate by LiH lines Al^h is given by 

_ Auh(LTE) 

- l + n„/n(H)- 
The critical density for LTE is 

Alih(LTE) 

"°^"'"'° AL.H|n(H)-.0| - (^°) 

where ALiH(LTE) is the cooling rate at the LTE, ALini'T'lH) 0] is that in the low density limit. 
For AuH[n(R) ^ 0], we use the fit by Galliand Palla (1998). 

2.1.3. Chemistry 

In addition to the hydrogen chemical network in paper I, we solve major deuterium 
chemistry between D, D+, D~, HD, and HD"*", and lithium chemistry between Li, Li"*", Li~, 
LiH, LiH"*". The included D reactions and their coefficients are the same as in Nakamura and 
Umemura (2002; Table 1). The Li reactions include those in Stancil, Lepp, and Dalgarno (1996) 
and are given in Table 1. 

Stancil, Lepp, and Dalgarno (1996) discussed possible enhancement of LiH via the three- 
body reactions 

Li + H + H-^LiH + H, (11) 
Li + H + H2->LiH + H2, (12) 

at high enough densities. We also include these reactions. Without accurate measuremnt of 
the rate coefficients for these reactions, we estimate them by assuming the same reaction cross 
sections as those of the corresponding hydrogen three-body reactions (Palla et al. 1983) but 
with correction owing to different molecular weight. The dissociative reaction 

LiH + H->Li + H2 (13) 
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Table 2. Initial model parameters. 



model 


T(K) 


y(H2) 


SI (fiducial) 


200 


5 X IQ-^ 


S2 


100 


5 X 10"'' 


S3 


300 


5 X lO""' 


S4 


400 


5 X 10-"* 


S5 


200 


2 X 10-3 


S6 


200 


10-3 


S7 


200 


10-4 



Table 2. The initial condition is taken at the number density lO'^cm ^ 

is important for limiting LiH fraction owing to its high reaction coefficient /cdiss = 2.0 x lO^^^cm^. 
We include the inverse reaction 

Li + H2 ^ LiH + H (14) 

with its rate coefficient calculated from the detailed balance. 

2.1.4- Luminosity Calculation 

For the thermal evolution, we employ an assumption similar to that for the dynamics: 
physical quantities evolve only when the mass element is inside the central region. After the 
mass element is left behind the central collapse at some instance, the physical quantities of the 
matter are kept constant until the end of the collapse phase. In this way, we can construct the 
distribution of temperature and chemical compositions in the entire cloud, as well as the density 
and velocity structures. Using this cloud structure, we calculate the molecular line emissivity 
in each part of the cloud (see 2.1.2) and estimate the total line luminosity by summing those 
from the central region and from the envelope. 

2.2. Initial conditions 

We start calculations at the number density of lO^cm"^. We calculate seven models 
whose initial temperature and H2 concentrations are given in Table 2. The model SI, which 
we call the fiducial model hereafter, represents the typical parameter set of first-star forming 
cores (e.g., Bromm et al. 2002). Around the fiducial parameters, we put possible variations of 
initial temperature and H2 abundance for other models. 

In particular, the difference in H2 abundance refiects the past history of star-forming 
cores in the parent cloud. If mass of the parent cloud is ~ IO^Mq, i.e., at the size of the 
first luminous objects in the standard ACDM cosmology, the H2 abundance reaches at most 
y{R2) ~ 10-^ — 10-'^ at the formation of star-forming core (~ lO^cm"^; e.g., Abel et al. 2002; 
Bromm et al. 2002). On the other hand, if the parent cloud is a giant pre-galactic cloud, 
i.e., ^ IO^Mq, and its virial temperature is higher than lO^K, the abundant H2 is formed 
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{y{B.2) ~ 10~^ — 10~^) owing to shock ionization and subsequent non-equilibrium H2 formation 
(Shapiro, Kang 1987; Susa et al. 1998; Oh, Haiman 2002). Considering the difference in 
parent cloud mass, we calculate models of both low- (models SI, S7) and high- (models S5, S6) 
H2 abundance. Similarly, for observing dependence on different initial temperatures from the 
fiducial one, we calculate models S2-4. The ionization degree is set to be y{e~) = 10~^. The 
abundances of other chemical species are fixed to the values calculated for the study of molecule 
formation in the primordial gas by Galli and Palla (1998): y(D) = 4 x 10"^, y{BD) = 1.2 x 10~^, 
y{U) = 10-^°, y(Li+) = 1.2 x 10"^°, y{UE) = 7.1 x lO^^o^ y{UE+) = 9.4 x 10"^^ The initial 
abundance of H^, Hf, H~, D"^, D~, HD"*", Li~ are set to zero. The concentration of the ith 
species y{i) is defined by y{i) = n{i)/ni{, where n{i) is the number density of the ith species 
and TiH is the number density of hydrogen nuclei. 

2.3. Numerical Results 

2.3.1. Evolution in the Fiducial Case 

The evolution of the star- forming core in the fiducial case (model SI) is presented in 
Figure 1, where (a) the temperature and (b) the H2, HD, and LiH fractions are shown against 
the number density. Here, the molecular fraction denotes the fraction of each element converted 
into the molecular species. For example, for LiH, the molecular fraction /(LiH) is defined as 
/(LiH) = n(LiH)/nLi, where nu is the number density of Li nuclei. In Figure 1(a), the solid 
curves show the evolution of the central region in the collapse phase. For the accretion phase, 
the evolution of three representative mass-shells is indicated by long- and short-dashed and 
dotted lines. The bold and thin lines indicate the temperature of the fiducial model and of the 
model with same parameters but without HD cooling, respectively. In Figure 1(b), the H2, HD 
and LiH fractions at the central region in the collapse phase are shown by solid and dashed 
and dotted lines. 

The small enhancement of HD fraction and resultant HD cooling in rK lO^cm"'^ (Figure 
lb) causes a slight deviation in temperature evolution from the case without HD (see Figure 
la). With increasing temperature for n > lO^cm"^, this HD is soon dissociated. The reason 
for this sudden HD dissociation is as follows: the abundance of HD in this density range is 
determined by balance between formation via 

D+ + H2^H+ + HD, (15) 

and dissociation via 

HD + H+^H2 + D+. (16) 

The equilibrium abundance between these reactions is given by y{B.D)/y{B.2) = 
2exp(464K/T)|/(D"'")/2/(H"''), which means that for low temperatures (T<464K), HD fraction is 
enhanced relative to H2 fraction. When the density exceeds the critical density ricr ~ lO^cm"^ 
for HD to reach the local thermodynamic equilibrium, the temperature begins to increase. 
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Since the increase in temperature reduces the HD fraction, HD coohng becomes inefficient. At 
higher density, the two evolutionary paths in Figure 1(a) converge. Accordingly, the presence 
of HD does not alter the evolution in the fiducial case for higher densities. In particular, the 
HD cooling does not modify the temperature evolution of mass shells in the accretion phase 
(Figure la). 

Hydrogen becomes fully molecular at ~ lO^^cm"^ via the three-body reactions. 



which start at ~ lO'^cm (Palla et al. 1983). Simultaneously, all the deuterium becomes HD 
via the molecular exchange reaction 



On the other hand. Figure 2 shows the effect of different LiH formation rate coefficients on the 
LiH fraction. Along with the case with the standard rate given in Table 1, following cases are 
shown: 

[Case 1] the molecular exchange reaction (reaction 1 in table 1) is excluded from chem- 
istry. The rate of three-body reactions (reactions 2,3 in Table 1) are the same as the standard 
values (case la) or is increased by a factor of ten (lb). 

[Case 2] the rates of the molecular exchange reaction and its inverse are reduced by a 
factor of 10 (case 2a) or 100 (2b) from the standard ones. In Figure 2, the LiH fraction in 
n < lO^^cm"'^ increases if the rate coeffecients of the reaction 1 and its inverse are decreased 
(i.e, case 2): since the LiH formation is suppressed by the dissociative reaction (equation 13; 
Stancil et al. 1996; Bougleux, Galli 1997; Puy, Signore 1998), lowering the rate coeffecient 
for this reaction results in enhancement of the LiH fraction. However, even in these cases, 
the LiH fraction remains < 10~^ for n < lO^^cm"'^ and their cooling rate is always negligible. 
At ~ lO^^cm"'^, the LiH fraction increases with the onset of the molecular exchange reaction 
(reaction 1 in Table 1). The Li three-body reactions are negligible in comparison with the 
exchange reaction for LiH formation. For higher density, the molecular fraction of lithium 
saturates at the equilibrium value ~ 10"'^ (Figure lb, 2). As in Figure 3, the total cooling rate 
by LiH L^n is less than lO^^erg s~^ in n < lO^^cm"^. At n ^ lO^^cm"'^, the maximum value 
Lun ~ lO'^^erg s~^ is reached, which is lower than and Lhd by 2-3 orders of magnitude. 
Some authors (Lepp, ShuU 1984; Stancil et al. 1996) considered that in some late stages 
of the collapse all the lithium might be converted to LiH by the three-body reactions, and 
thus LiH could be an important coolant. However, in our calculation, LiH never become an 
important coolant in any evolutionary stage. As for HD, Lhd is smaller by more than an order 
of magnitude than for n > lO^cm"'^, although Lhd is comparable to Lhj at < lO^cm"^ 
in the fiducial case (Figure 3). 



3H^H2 + H, 
2H -|- H2 — > 2H2, 



(17) 
(18) 



D + H2 



H + HD. 



(19) 
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2.3.2. Dependence on Initial Conditions 

In the following, we describe the dependence on initial conditions. First, we show the 
dependence on initial temperature by comparing the fiducial model with models S2-4 (Figure 
4). For the initial temperature ^ 200K, temperature drops, since the H2 cooling is effective in 
this temperature range. On the other hand, for lower initial temperatures, the temperature rises 
with the compressional heating owing to inefficient cooling. Then, the temperature approaches 
around T ~ 200K, and for n ~ lO^cm"^, these evolutionary trajectories converge in all these 
cases. 

Next, we see the dependence on initial H2 abundance yi(H2) by comparing between the 
fiducial model and models S5-7 (Figure 5). For low initial H2 abundances yi{^2) ^ 10~^, the 
temperature does not drop sufficiently for HD formation. In this case, without enough HD 
formation, the temperature evolution follows the fiducial case closely (see Figure 5a). The 
situation changes for higher values of initial H2 abundance ?/i(H2) > 10~^, (Uehara, Inutsuka 
2000; Nakamura, Umemura 2002). With such a high value of ?/i(H2), the temperature falls 
below the threshold value (~ 150K) for HD formation. Because of the rapid HD formation 
and consequent HD cooling, the temperature drops further. Although HD cooling becomes 
inefficient after the density exceeds the critical value for LTE (~ lO^cm""^), the temperature 
evolution differs from the fiducial case owing to different H2 abundance. When the difference 
in H2 abundance becomes negligible due to the three-body reactions at ~ lO^cm"^, the tem- 
perature evolution converges. In above cases with high initial H2 abundance (models S5 and 
S6), the temperature rises adiabatically around n = 10^ — lO'^cm"'^. In this regime, the actual 
contraction of clouds should be slower than the assumed free-fall rate. Longer contraction 
timescale reduces the compressional heating rate, and thus the temperature increase should be 
more gradual than in our model. Nevertheless, alteration of the evolution at this density range 
does not affect the maximum luminosity: in such high density regime where the maximum 
luminosity is attained, the evolution has converged already regardless of the behavior in the 
lower density regime. 

In conclusion, the evolutionary trajectories of both temperature and HD abundance are 
most affected by the initial y(ii2) value. This difference disappears n ~ lO^cm"'^ with rapid H2 
formation by the three-body reactions. The evolution at high densities is independent of the 
initial conditions. In particular, the evolution of inner mass shells {Mr ^ SOMq) is independent 
of the initial conditions because the density there is higher than lO^cm"^. 

2.3.3. Luminosity Evolution 

The maximum luminosities Lpcak of the two strongest H2 lines 1-0 Q(l) and 0-0 S(3) 
^, and the mass of the protostar when these maximum luminosities are attained, Mpeak are 

^ For example, the notation 1-0 Q(l) means the following transition: " 1-0 " denotes vibrational transition 
from u = 1 to 0. Q denotes the change of rotational quantum number AJ — 0. Similarly, O, P, R, and S, 
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presented in Table 3 for each model. 



Table 3. Maximum Line Luminosities. 



model 


i-noakfl-O QfDldO^^ere; s"^) 

pcaK V ^ V / / V o / 


'-^ pcaK \ \z) J 


L„™k(0-0 S(3))('103W s~M 




SI 


1.69 


12.9 


1.11 


18.2 


S2 


1.70 


13.8 


1.11 


18.2 


S3 


1.69 


13.8 


1.11 


18.2 


S4 


1.69 


13.8 


1.11 


18.2 


S5 


1.50 


12.0 


0.70 


7.9 


S6 


1.63 


12.9 


0.91 


10.5 


S7 


1.94 


29.5 


1.33 


19.5 



As examples, we choose the two strongest H2 lines in our calculation, the rovibrational 
line of 1-0 Q(l) and pure rotational line of 0-0 S(3), whose wavelengths at rest are 2.34/im and 
10.03/im, respectively. For models SI, S2, S3, and S4 (cases with same initial H2 abundances 
but different initial temperatures), the luminosities evolve in a similar way, since both the 
temperature and HD abundance converge immediately. In Figure 6(a), we show the luminosity 
evolution of models SI, S6, and S7, which are cases with the same initial temperatures, but with 
different initial H2 abundances. The temperature, H2 and HD abundances in these models are 
significantly different up to lO^cm"^. For lower densities, the luminosity from the cores depends 
on initial conditions (Figure 6a), while for higher densities (> lO^cm"^) physical quantities 
converge. Since temperature is higher in this range than in lower densities and the emission 
from the highest temperature region dominates the total emission (see §4.1 later), the similar 
temperatures lead to convergence of H2 line luminosities for both the late collapse phase and 
the early accretion phase. Note that regardless of initial conditions, the peak values of the 
rovibrational line luminosities are almost unchanged: among our models, variations in |/i(H2) 
cause difference in the peak values by less than a factor of two. In Figure 6(b), we show the 
luminosity evolution of the HD lines 0-0 R(l) (56.2/im) and 0-0 R(7) (15.0/xm), and that of 
the H2 lines 1-0 Q(l) and 0-0S(3) for the model S6. This is the case with the same parameters 
as in paper I. Although the main coolant is still H2 as in paper I, its cooling rate is reduced 
because of the contribution by HD cooling for the early collapse phase. Since physical quantities 
converge for higher densities (> lO^cm"^) regardless of HD contribution in low density region 
(Figure 1, 5), the maximum luminosities of H2 lines (Table 1 in paper I) remain almost the 
same. 

denote A J = 2,1,-1, and -2, respectively. The number in brackets denotes the rotational quantum number 
after the transition. 
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3. Evolution of a Filamentary Cloud 



A large pre-galactic cloud IO^Mq) collapses into a sheet-like configuration because of 
an instability against non-spherical perturbation (Lin et al. 1965; Hutchins 1976; Susa et al. 
1996). This oblate cloud experiences a shock and forms a postshock-cooled layer. The cooled 
layer more likely fragments into filamentary clumps than into spherical ones (Miyama et al. 
1987a, 1987b). In the shocked layer, formation of H2 and HD is promoted because of abundant 
free electrons produced by the shock heating, which catalyze the formation reactions (MacLow, 
ShuU 1986; Shapiro, Kang 1987). The filaments continue to contract due to effective molecular 
cooling. At ~ lO^cm"^, where heat input by H2 three-body formation becomes significant, the 
filaments fragment and dense star-forming cores are produced (Nakamura, Umemura 2002). 
This is in contrast to the situation in a small cloud (~ IO^Mq), where cores directly form at 
~ lO^cm-3 (Bromm et al. 1999,2002; Abel et al. 2000,2002). Because of a larger molecu- 
lar abundance, stronger molecular emission is expected in filamentary clouds. Here, we then 
evaluate the molecular line luminosity from those filaments. 

3.1. Model 

In studying the thermal and chemical evolution of a filamentary cloud, we approximate 
the filament by an infinitely long cylindrical cloud, for simplicity. Since a sheet fragments when 
pressure retards gravitational contraction, in filaments thus formed, the pressure is comparable 
to the gravity. Moreover, the ratio of gravity to pressure force is constant during the isothermal 
evolution, which is indeed approximately the case. The small deviation from the isothermal 
evolution has an important dynamical effect. Here, we need to treat the effect of gas pressure 
carefully in our model. 

A contracting filament consists of two parts, a denser spindle with almost flat density at 
the center and an envelope with lower denisty (e.g., Inutsuka, Miyama 1997). As in the case of 
spherical cores (section 2), we focus on this central dense spindle, and apply the virial equation 
for a cylinder with uniform density: 



2 dt^ 



2T + 2Yi-Gm\ (20) 



where 



X= I pr^dV, T= f -pu'^dV, U= f pdV, (21) 
Jv Jv 2 Jv 

where the integration acts over the volume per unit cylinder length V , and m is the line density 
of the central spindle. Assuming constant density p in the spindle and using scale radius 



R = ^Jm/np^ the integration can be worked out: 

1 1 .dR.r, „ kfiT , . 

Z = -mR^, T = -m(—?, H = -^m, 22 
2 ' 4 ^ dt^ pmn 

where p is the mean molecular weight, and mn is the mass of a hydrogen nucleus. Substituting 
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equation (22) into equation (20), we obtain 
2G 



m - m^.itiT)], (23) 



c/t2 R 
where the critical hne density 

Note that the gas pressure exactly balances the gravity for the critical line density merit (7"). 

The line density of the spindle m is not constant during the collapse. If the gravity is 

far larger than pressure, the outer part of the spindle is left to the envelope. In other words, 

when the line density of the central spindle m is much larger than the critical value merit, ^ 

decreases during the collapse of the cylinder. To mimic this effect, we assume that if m is larger 

than a critical value, m decreases at the same rate as the core radius does: 

dm 1 dR , , 

^ = (m-mtir)^^. (25) 

Here, we take the threshold line density mth = l.Smcrit following Nishi (2002), who found this 
treatment reproduces the numerical results of Nakamura and Umemura (1999) well. 

As for the envelope, we employ an assumption similar to the case of spherical cores 
(section 2.1.1): physical quantities of a mass element evolve only when the mass element is 
inside the central spindle, and remain constant after the mass element is left to the envelope. 
Therefore, by calculating only physical quantities of the central spindle, we can obtain the 
whole cloud structure. 

Thermal and chemical processes are treated in the same way as in the spherical case. 
The luminosity per unit length is also calculated in a way similar to before, by summing the 
emissivity both in the dense spindle and the envelope. When we compare this with the spherical 
ones, we set the filament mass to be Mfn = IO^Mq, and thus length, /gi = Mfn/nii. 

When the fragmentation timescale of filaments tfrag = 2.1/ {2iiGpY^'^ (Nagasawa 1987) 
becomes shorter than the collapse timescale tdyn = p/{dp/dt), fragmentation is expected to 
occur. We follow the evolution of the filaments up to this epoch. 

3.2. Initial condition 

Collapse and fragmentation of a filamentary primordial gas cloud have been investigated 
by some authors (e.g., Uehara et al. 1996; Uehara, Inutsuka 2000; Nakamura, Umemura 2002). 
Following the classification of the filaments by Nakamura and Umemura (2002), we study three 
varieties of filaments: 

• Fl: a low-density filament with low H2 fraction, whose contraction is controlled by II2- 
cooling: the initial number density = lO^cm""^, temperature T = 400K, and H2 abun- 
dance ?/(H2) = 10~^, 

• F2: a low-density filament with high H2 fraction, whose contraction is controlled by HD 
coohng: un = lO^cm'^, T = 300K, and y(R2) = 3 x 10"^ 
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• F3: a high-density filament: tt-h = lO'^cm"^, T = 400K, and y(H2) = 10"^. 

We use the initial ionization degree in filaments y(H^) = ?/(e"") = 5 x 10~^ following Nakamura 
and Umemura (2002). Except H2, H+, and e~, we use the same initial abundances of chemical 
species as in the case of the spherical cores. 

Using the ratio of gravity to pressure gradient /, and the critical line mass merit (^init) for 
initial temperature Tinit, the initial line density mi can be written as minit = fmcvitiTmit)- Here, 
we assume / = 2, which is the most probable value for filaments formed by fragmentation of a 
sheet-like cloud (Miyama et al. 1987a,b). For the total mass of lO^M©, the filament lengths 
are (Fl)32.5 pc, (F2)43.4 pc, and (F3) 32.5 pc. 

3. 3. Results 

3.3.1. Thermal and Chemical Evolution 

Figures 7-9 show evolution of (a) temperature; (b) abundances of H2 (solid line) and HD 
(dashed line); (c) cooling rates by H2 (solid line) and by HD Lhd (dashed line). We do not 
present the evolution of LiH abundance here, since LiH formation proceeds in a fashion similar 
to the case of the spherical cores: only a tiny fraction of Li becomes molecular form without 
affecting thermal evolution. 

The overview of evolution in each case is as follows: 

• Fl: the low-density filament with low H2 fraction 

The evolution resembles that of the spherical cores with low H2 fraction (e.g., model SI in 
section 2). HD fraction remains low (Figure 7b). Throughout the contraction, dominant 
cooling mechanism is H2 line cooling, which nearly balances the compressional heating 
(Figure 7c). In densities higher than the critical density for collisional de-exitation of H2 
rotational level (~ lO^K), temperature rises slowly from 400K to 700K (Figure 7a, c). With 
the onset of rapid H2 formation via the three-body reactions and resultant heat injection, 
the temperature rises suddenly and the pressure becomes stronger than the gravity (Figure 
7c). At this moment (~ 10^'^cm^'^) the fragmentation timescale tfrag becomes shorter than 
the collapse timescale tdyn and then the cylinder fragments. 

• F2: the low-density filament with high H2 fraction 

As in the case of the star-forming cores with high H2 fraction (model S5, 6 in section 2), 
HD is formed abundantly and becomes an important coolant (Figure 8a,c). As a result, 
the temperature goes down to < lOOK. When the filament contracts n ^ lO^cm""^, which 
is the critical density for HD collisional de-excitation, the temperature begins rising as 
described in section 2.3.1. Owing to this temperature increase, HD begins to dissociate. 
This results in further temperature increase and the fragmentation occurs at < lO^cm"'^. 

• F3: the high-density filament 

The initial density being already higher than the LTE critical density for H2 and HD 
(~ lO^cm"^ for H2 and ~ lO^cm"^ for HD), the temperature increases with increasing 
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density. 

This filament begin to collapse at rather high density from a highly gravitationally unstable 
(/ = 2) state. Even in a density range as high as ~ lO^cm"^, the gravity dominates over 
the pressure, and the filament contracts faster than in other models. Then, in spite of H2 
formation heating around n ^ 10®cm~^, this filament contracts without fragmenting until 
~ lO^^cm"'^ where it becomes opaque to the H2 lines. 

3.3.2. Luminosities from the Filaments 

In Figure 10, time- averaged luminosities of H2 and HD lines from the filamentary clouds 
are shown. Here, time averaging is over the total contraction time before fragmentation. Total 
mass of the filaments is normalized to IQ^Mq. 

In model F3 (Figure 10c), the filament evolves in the region of high temperature 
(400K^ T ^1500K) and high density (lO^cm"'^ ^n'^ lO^^cm"'^) where H2 cooling is efficient. 
Accordingly, strong radiation is emitted with some H2 lines exceeding lO'^^erg s~^ and HD lines 
exceeding 10'^^ erg s~^. Luminosities are weaker for model Fl where all the line luminosities 
are less than lO'^^erg s~^ (figure 10a) because of the lower temperature (400- TOOK) in model 
Fl throughout the contraction. With temperature even lower (T <200K), in model F2 only a 
few lower-lying rotational lines of H2 and HD are stronger than lO'^^erg s~^. Since HD is more 
effective as a coolant than H2, most of the H2 lines are very faint (< lO^^erg s~^). 

For the purpose of seeing the effect of cloud geometry, we show in Figure 10(d) the mean 
luminosities of spherical star-forming cores of total mass lO^M© calculated with the same initial 
condition as for model (F3). We can see that line luminosities of model (F3) are stronger than 
those of star-forming cores with the same mass by about an order of magnitude. The reason 
is as follows: In the case of a spherical core, the fully molecular region, where H2 emission is 
strongest owing to heating by H2 formation, accounts for at most ~ IMq at the center (Omukai, 
Nishi 1998). Most of the matter is left to the envelope before this epoch. On the other hand, in 
the cylindrical case, the fraction of envelope mass is less than that in the case of the run-away 
collapse phase. Although the radius of the fully molecular region at the center of the cylinder 
is almost the same as in the spherical case, its length in the axis direction is longer. Therefore, 
the mass of the fully molecular spindle is larger and emission from it is stronger. 

4. Discussion 

4.1. Analytical Expression for the Luminosity 

In this section, we derive an analytical expression for the luminosity from the collapsing 
clumps, which reproduces well our numerical results. First, let us consider the luminosity from 
a spherical core in the collapse phase. During the run-away collapse phase, length scale of the 
central region is about a Jeans length Aj and mass of the region M^en is about a Jeans mass 
Mj = pAj. Using cooling rate (per unit volume) in the central region Accn, the luminosity from 
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the central region is given by 

Lccn — -^ccn-Accn/ Pcen- (2^) 

Here we assume an energy balance between heating and cooling, which is valid in a good 
accuracy during the run-away collapse phase: 

■^ccn — Teen- (^'^) 

Also, we consider only the compression as a source of heating, then 

/xmn iff 

where we have assumed that the timescale of contraction is approximated by the free-fall time. 
Substituting equations (27) (28) into (26), the luminosity from the central region is given by 

/ Q2 / T \ ^ 
Lcen^ J — 7^^46L0 ergs"^ (29) 

Note that L^cn is a function only of temperature: it does not depend, for example, on mass or 
density of the central region. In other words, if the temperature is constant throughout the 
evolution, the luminosity from the central region remains constant. 

To calculate the total luminosity, we need to add the contribution from the envelope. 
In the collapse phase, the Jeans mass decreases with increasing density. A mass shell initially 
contained within the central region at some time leaves the central evolution and becomes part 
of the envelope when the mass inside the shell equals the central Jeans mass. We assume that 
physical quantities in the envelope are frozen: physical quantities on a mass shell are kept 
constant after the shell leaves the central region until the end of the collapse phase. Then, the 
cooling rate of a mass shell in the envelope is given by that of the central region at the moment 
when the shell is detached from there. Because of homogeneity inside the central region, the 
cooling rate in this region is 

(-] = #^ (30) 

where L^cn and M^cn is the luminosity and mass of the central region. Since a shell that contains 
mass M inside is detached from the central region when M^en = M, the cooling rate of the mass 
shell in the envelope A(M) is given by 

A(M) = ^^. (31) 

Here we use the constancy of luminosity of the central region i.e., Lccn{M) = L^en (see equation 
29), considering a star-forming core collapsing isothermally. Summing over all the mass shells 
in the envelope, the luminosity of the envelope is 

/•Mtot A 

Lenv=/ -{M)dM (32) 
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.'"(■^). (33) 



cen 



where Mtot is the total mass of the star-forming core. Then, the total luminosity of the star- 
forming core is 



(34) 

= LceJl + \n(^)]. (35) 

The envelope contribution enhances the total luminosity only logarithmically: with increasing 
central density and decreasing M^cn, total luminosity increases only very slowly. In the following, 
we do not discriminate between Lccn and Ltot- Recall that, in deriving the expression above, we 
have not specified a mechanism for cooling. In the case of primordial gas, which we consider 
here, cooling is via H2 and HD line cooling. If many lines are equally important for cooling, the 
total luminosity is divided into them. However, in many situations, one or two lines become 
dominant. In such cases, the luminosities of these lines are given approximately by the total 
luminosity Lt^t- From equation (29), we also see that when HD becomes more important than 
H2 in cooling, low temperature (~ lOOK) causes a corresponding drop in luminosity. 

Here we compare the analytical expression (29) with the numerical results in section 
2. In Figure 11, we show the H2-line luminosity for model SI and the analytical values 
(29) for the model SI temperature evolution. Also shown are the analytical values for constant 
temperatures T = 500K and lOOOK. For n < lO'^cm"^, the numerical value matches well 
the analytical value Lccn- On the other hand, the numerical result deviates from the analytical 
value for higher densities {n > lO^cm^'^) because the assumption that the main heating source 
is gravitational compression, which is postulated in deriving equation (29), is violated: for 
these high densities, dominant heating process is the chemical energy released by H2 formation. 
In particular, the strongest emission comes from the region of n = 10^ — 10^°cm~^, where H2 
formation is most active. For higher densities {n > 10^°cm~^), since majority of hydrogen has 
already become the molecular form, the chemical heating decreases and has little contribution 
to the entire luminosity. Then, when the central density is larger than 10^°cm~^, the total 
luminosity is dominated by the contribution from the envelope of n = 10^ — 10^°cm~^. The 
constancy of luminosity for n > 10^°cm~^ is the result of the dominance of the emission from 
the envelope, which hardly evolve during the collapse of the center. 

Next, we discuss the case of filamentary clouds. The size of the central spindle is again 
about a Jeans length Aj and the line density of the spindle is ~ Pccn-^j- The contribution of 
luminosity from the central spindle is then 

-^^sp - PccnAj/(Accn/Pccn), (36) 

where / is the filament length. As in the case of spherical cores, the thermal balance (equation 
27) is maintained during the collapse. If heating is due only to compression, using equations 
(27), (28) and (36), the luminosity of the spindle is given by 
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- P«..^lc J^i^ ^ 4.2 X lO'Le (^)' (^) * (^) . (37) 

Unlike spherical cores, the luminosity of a filamentary cloud rises with increasing density even 
in the isothermal case. So far, we have considered only the contribution from the central 
spindle. In the isothermal case, since the line density of the central spindle, pccnAj, is constant, 
all the matter in the filament contracts uniformly without being left to the envelope. The 
central contribution Lgp becomes equal to the total luminosity. Indeed, the filaments remain 
at approximately constant temperature in many cases in section 3. Then, we assume in the 
following that Ltot — Lm. 

Next, we compare the analytical expression (37) with the numerical results. We consider 
the range n < lO^cm"^, where H2 formation heating is not important. In figure 12, luminosity 
evolution of models Fl-3 is shown, along with analytical values shown by equation (37) for 
T = lOOK, 500K, and lOOOK. The temperatures in models Fl and F3 are about 500K and 
remain roughly constant, while in model F2 the temperature is about lOOK, but with significant 
variation. For model Fl, the analytical value matches the numerical one by a factor < 3. 
Agreement between them is even better for model F3. In the case of model F2, the numerical 
luminosity deviates from the analytical one for lOOK because of a large variation in temperature. 
When we calculate the analytical value using the temperature evolution of the numerical result, 
our numerical result agrees with the analytical one within a factor of a few. In summary, the 
analytical estimates reproduce our numerical calculation within a factor of a few. 

4-2. Expected Luminosity from Pop III Star- Forming Galaxies 

In this section, we discuss the condition of the Pop III star-formation rate (SFR) ipm 
and formation redshift z for enabling the detection of molecular emission by next generation 
facilities, in particular, JWST (optical-mid-infrared wavelengths) and SPICA (mid-far-infrared 
wavelengths). Hereafter, we consider Pop III star-forming galaxies where star formation pro- 
ceeds through low-density filaments or through high-density filaments. 

First, we discuss the case of star formation through the low-density filaments. In this 
scenario, in a large pre-galactic cloud, shocked layers formed at virialization fragment into low- 
density (~ lO^cm"^) filaments, which again fragment at ~ lO^cm"^ into star-forming cores of 
10 — 4OM0 (Nakamura, Umemura 2002). In the following, the mass of the cores is assumed to 
be ~ IOM0. Under the assumption that the relative number of cores in a specific evolutionary 
phase is proportional to the duration of that phase, the luminosity averaged over the lifetime 
of a core Lcore coincides with the ensemble mean luminosity of cores in a pre-galactic cloud. 
Here, the physical quantities at the fragmentation of the low density filaments are similar to 
those of the star- forming core at ~ lO^cm"'^ in model S6 of section 2. Then, in the following 
numerical evaluation, we use the model S6 core at lO^cm"^ as a model of the fragments. 
The duration of the core collapse is ~ Atcore = lO^year and the time-averaged luminosities of 
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the strongest lines from the cores, which include H2 rovibrational lines of 1-0 Q(l)(2.34/im), 
1-0 0(3)(2.69/im), 1-0 O(5)(3.04/im), and H2 pure rotational lines of 0-0 S(3)(10.03/im), 0-0 
S(4)(8.27/im), 0-0 S(5)(7.07/im)), are Lcorc ^ (2 - 6) x lO^^erg s^^ On the other hand, for 
the low-density filament before fragmentation into cores (model F2 in section 3), the strongest 
lines are H2 pure rotational lines: 0-0 S(0)(29.7/im), 0-0 S(l)(17.9/im), and HD pure rotational 
lines: 0-0 R(0)(112/im), 0-0 R(l)(56.2/xm) with Zm ~ (2-5) x lO^^erg s~^ The duration of the 
filamentary contraction is roughly Atcorc = 4 x lO^year. Since the strongest lines differ between 
the core and the filament, we can judge the conditions of emission regions by the line diagnosis. 
Next, we estimate luminosity from the whole pre-galactic cloud using the values of luminosity 
from the individual cores/filaments given above. With luminosity of a line from the core Lcore 
and number of cores in cloud Ncom, the luminosity of the line from the whole cloud Ld^com is 
described by Ld^core = -^coro^core- The total number of star formation cores A^coro is ipiw^t/M^, 
using the Pop III SFR ipm, the typical mass of formed star M^, and the lifetime of the core At. 
Here we assume that the star formation efficiency e is unity because the density is sufficiently 
high (n ~ lO^cm"'^) for star formation. Using the luminosity of the strongest lines (e.g., H2 
0-0S(3) at 10.03yum or 0-0S(5) at7.07/im) from an individual core Lcore = 5 x 10^^(erg s~^), the 
stellar mass M^, = IOMq, and the lifetime of the core At = lO^yr, the luminosity from the whole 
pre-galactic cloud is 



Next we consider the contribution from filaments. The number of filaments in a pre-galactic 
cloud is A^fii = -j/^/j/Atfii/eMfii. In numerical evaluation, we use e = 0.1, which is a typical value 
of star formation efficiency of some cluster forming regions (e.g., Megeath et al. 1996; Kumar 
et al. 2004). Using the contraction time of filaments Atai = 4 x lO^yr, the mass of a single 
filament Mai = IO^Mq, the luminosity of the strongest line (i.e., HD O-OR(O) at 112/im) of a 
single filament Lgi = 5 x lO'^'^erg s~^, the sum of luminosities from filaments is 



Thus, the flux from the cloud at redshift z 



where is the luminosity from the whole pre-galactic cloud and -D1+2 is luminosity distance at 
redshift z. We use the cosmological parameters Hq = 70kms~^Mpc^^, Q\ = 0.7 and Qu = 0.3. In 
Figure 13 we show the redshift-SFR domain where the H2-line flux exceeds the detection limit 
of SPICA, 10~^^ — 10~^^Wm~^ in the mid- and far-infrared wavelengths with a spectroscopic 
resolution R = 3000 and observation time of 10 hours (private communication; H. Matsuhara). 
SPICA is the most sensitive in the redshifted wavelengths of molecular emission for 2; < 10 among 
next-generation facilities. In Figure 13, the top and bottom lines are drawn by the criterion 




(38) 




(39) 
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F = Lci/4:'kDI_^_^ > 10~^^Wm~^ using the luminosities from the whole cloud Ld = lO^'^^ipmeig s~^ 
and = 5 X lO^'^i/jmeig s~^, respectively. With Pop III SFR as high as ~ 10^ — lO^Moyr"^ at 
redshift 2; < 8 — 10, SPICA can detect the strongest H2 and HD lines, such as H2: 1-0 Q(l), 
1-0 0(3), 1-0 0(5), 0-0 S(0), 0-0 S(l), 0-0 S(3), 0-0 S(4), 0-0 S(5), HD: 0-0 R(l). Being 
more sensitive than SPICA for A < lO/im, JWST can detect H2 vibrational lines under the same 
condition in Figure 13. However, strong H2 and HD pure rotational lines are redshifted away 
from the JWST band for 2; ^ 3. To study the physical condition of emission regions using the 
luminosity ratio of rotational and vibrational lines, SPICA is more suitable. 

Next, we discuss the case of star formation through the high-density filaments (^ 
lO^cm"^). Because the formation mechanism of these high-density filaments remains unclear, 
we do not try to include the emission from the lower density region and consider only the con- 
tribution from the high-density filaments (model F3 in section 3). The strongest emission lines 
in this case are very similar to those from spherical cores, and their time-averaged luminosities 
are Lgi ^ lO^^erg s~^ for filament mass of IO^Mq (see figure 10c). The number of filaments in 
the pre-galactic cloud Nm = (^/^///Atfii/eMfii). Because all the matter in the filaments is dense 
enough to form stars, we assume that SFR in these filaments is unity: e = 1. The sum of 
luminosities in a strong line over all the filaments in the pre-galactic cloud 



where we have used the contraction time Atgi = lO^yr and the mass of a filament = IO^Mq 
and the luminosity of the strongest lines (e.g., H2: 1-0 Q(l), 1-0 0(3), 1-0 0(5), 0-0 S(3), 0-0 
S(4), 0-0 S(5)) Lfii = lO'^^erg s~^. The z — tfjin domain for detection of the molecular emission 
by SPICA for the high-density filaments is also shown in Figure 13. We can see that SPICA and 
JWST can detect the molecular emission if the SFR is as high as ~ 10^~^MQyr~^ at redshift 



In summary, in both cases, high (~ lOOOMoyr"^) Pop HI SFR at redshift 2 < 7 is required 
for protogalaxies to be observable by next-generation facilities through the molecular emission. 

How much mass is needed in these pre-galactic clouds ? Duration of star formation 
in a pre-galactic cloud would be larger than the free-fall time of the cloud, tg ^ lO^yr for 
~ lcm~'^. For star formation to continue at a rate of ipm during this period, gas mass of 
ipiiits is required. For ipm > IO^MqJt'^, the gas mass of a pre-galactic cloud exceeds IO^^Mq. 
Also, metallicity in the cloud must be lower than 10~^ — 10~^Zq to cool and emit a copious 
amount of primordial molecular emission (Omukai et al. 2005). In the standard bottom-up 
scenario of galaxy formation, when building blocks assemble into a protogalaxy as massive as 
~ 1O^M0, the volume-averaged metallicity of gas has reached [Fe/H] = —2.30 via metal diffusion 
by supernovae (Mori et al. 2004). Their results also show that the metal mixing proceeds quite 
imhomogeneously. Large spatial variation of metallicity ranging from [Fe/H] — 5 to is found 
in high-density (> lcm~^) regions. While the ambient gas of supernovae is polluted by metals. 




(41) 



z<7. 
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most of the gas remains at low metallicity, [Fe/H] ~ —5. Within a relatively metal-polluted pre- 
galactic cloud, a large number of almost metal-free regions can be contained. If star formation 
is active in such very metal-poor regions, next-generation facilities might have a chance of 
detecting the H2 and HD lines emitted in those regions. 

4-3. Observational Strategy for Discovering Pop III Galaxies 

Finally, we propose a strategy for discovering pre-galactic clouds that contain Pop III 
star-forming regions. Such objects are likely to be metal- deficient. A method for finding such 
objects is then to look for the La-line emission, which is expected to be emitted copiously 
(~ lO^^erg s~^) at the virialization of massive (10^^ — IO^^Mq) pre-galactic clouds (Omukai, 
Kitayama 2003). A number of La blobs and La emitters have been found in recent surveys. 
Although the emission mechanism of those objects still remains disputed, among hypotheses 
is the cooling radiation at the virialization (Steidel et al. 2000, Haiman et al. 2000, Fardel et 
al. 2001). Some La blobs are very diffuse and in their vicinity no bright continuum source is 
found. These properties are indeed expected for the La-cooling objects (Matsuda et al. 2004). 

Recently, La blobs that are brighter than lO^^erg s~^ and larger than lOOkpc have been 
found in and around a protocluster region at redshift z = 3.1 (e.g., Steidel et al. 2000, Matsuda 
et al. 2004). Those La blobs meet the detection condition derived in section 4.2 above, namely, 
being Milky- Way sized and located at 2 ~ 3. If those La blobs indeed originate in the cooling 
radiation of low-metallicity gas, there might be a number of Pop III star-formation regions 
inside them. As an observational strategy, we propose first to select candidates by investigating 
metallicities of La blobs with spectroscopic analysis and then to look for the molecular lines 
using SPICA or JWST (see also Omukai, Kitayama 2003). 

5. Summary 

We have studied chemical evolution of primordial molecules (H2, HD and LiH) in star- 
forming clumps and evaluated their line luminosities. We have found conditions of the star 
formation rate and the formation redshift of a Pop III star-forming pre-galactic cloud to be 
observable by SPICA and JWST. 

In the case with high initial H2 abundance, HD cooling becomes effective at low density 
n < lO^cm"^ and affects the thermal evolution. At densities higher than the LTE critical density 
~ lO^cm"^, temperature increases and HD is dissociated. At high enough density (n > lO^cm"^), 
H2 dominates the cooling and HD has negligible effects on the evolution. LiH cooling remains 
insignificant despite the expectation that LiH can be an important coolant owing to abundant 
LiH formation at the high density {n > lO^^cm"^; e.g., Lepp, ShuU 1986; Stancil et al. 1996). 
Even though the three-body LiH formation is included, the molecular fraction of lithium remains 
low (< 10~^) owing to efficient dissociation through LiH -|- H ^ Li -|- H2. 

Molecular emission from a star-forming core reaches the maximum shortly after the 
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formation of a central protostar. Since the evolution has already converged by this moment, 
maximum luminosities of lines hardly depend on the initial conditions. Most luminous lines 
include H2 lines of 1-0 Q(l), 1-0 0(3), 0-0 S(3), and 0-0 S(4). Their luminosities are ^ 
lO^^erg s~^ for the core of IO^Mq . For a contracting high-density filament, strongest lines 
include H2 lines of 1-0 Q(l), 1-0 0(3), 1-0 0(5), 0-0 S(3), 0-0 S(4) and 0-0 S(5). The 
luminosities of these lines are ^ 10^''erg s~^ for a filament of 1O^M0. On the other hand, for 
a low-density filament of lO^M©, the luminosities reach > lO^^erg s~^ for the strongest lines 
( H2 lines of 0-0 S(0), 0-0 S(l), and HD lines of 0-0 R(0), 0-0 R(l) ) before fragmentation 
into cores. In protogalaxies with a high Pop III star formation rate ~ 10'^ — lO^MQyr"^ at 
redshift z<7, SPICA is able to detect H2 hues such as 1-0 Q(l), 1-0 0(3), 1-0 0(5), 0-0 
S(0), 0-0 S(l), 0-0 S(3), 0-0 S(4), and 0-0 S(5), and the HD line 0-0 R(l). Since JWST is 
more sensitive than SPICA in the near-infrared wavelengths ^ lO/im, JWST also can detect 
the H2 rovibrational lines. SPICA, whose bands include mid- and far-infrared wavelengths, has 
advantage of observing redshifted H2 and HD pure rotational lines emitted z>3. Then, we can 
reveal physical conditions of star-forming regions by diagnosing the rotational and vibrational 
line ratios. Provided such objects exist, these facilities can directly probe Pop HI star formation 
through H2 lines. 
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Tohru Nagao, Francesco Palla, Masayuki Umemura, and Torn Yamada for valuable comments, 
and Ken-ichi Oohara and Kazuya Watanabe for continuous encouragement. This work is sup- 
ported in part by research fellowships of the Japan Society for the Promotion of Science for 
Young Scientists and a Grant-in-Aid for Scientific Research (16204012 KO; 14540227 RN) of 
the Ministry of Education, Culture, Science and Technology. 
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Fig. 1. Evolution of (a) the temperature and (b) the molecular fractions of H (solid) and D (dashed) and 
Li (dotted) against the number density for the fiducial model (SI). The molecular fraction is the fraction 
of each element locked into the molecular species shown. In (a), the temperature in the central part in the 
run-away collapse phase is drawn with a thick solid line. Here, that of the model with the same parameters 
as in the fiducial model but without HD cooling (thin line) is also shown for comparison. For the accretion 
phase, the evolutionary trajectories of the mass shells oi M/Mq = 1 (long-dashed fines), 13 (short-dashed 
lines) and 50 (dotted lines) are shown. In (b), only the evolutionary trajectories of the central region are 
shown. The II2 fraction is shifted by 10^ for clarity. 




Fig. 2. Evolution of LiH abundance for different formation rate coefficients. Shown are the values at the 
center in the collapse phase. The standard case, whose rates are given in Table 1, is shown by the solid 
line. The cases la, b, where the molecular exchange reaction (reaction 1 in Table 1) are omitted, are shown 
by long-dashed (la) and short-dashed lines (lb). In the case la the coefficients of the three-body reactions 
(reactions 2,3 in Table 1) are the same as the standard values, while in the case lb they are increased by a 
factor of ten. The cases with reduced rates of the molecular exchange reaction and its inverse by a factor 
of 10 (case 2a) or 100 (2b) from the standard ones are shown by dotted and long-dash-dotted lines. 
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Fig. 3. Evolution of the cooling rate by H2 (solid line), HD (dashed line) and LiH (dotted line) in the 
run-away collapse phase against the number density for the fiducial model (SI). 
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Fig. 4. Effect of difference at the initial temperature Tinit: (a) temperature evolution and (b) evolution 
of the H2 (upper lines) and HD (bottom lines) fractions against the number density for the run-away 
collapse phase. The evolutionary trajectories of model SI (fiducial): Tinit = 200K (solid), S2: Tinit = lOOK 
(long-dashed), S3: Tinit = 300K (short-dashed) and S4: Tinit = 400K (dotted) are shown. 




4 6 8 10 12 14 4 6 8 10 12 14 



log[ n ( ] log[ n ( cm-3) ] 

Fig. 5. Effect of difference in the initial H2 concentration yinit(H2). Meaning of the panels is the same 
as in Figure 2. The evolutionary trajectories of model Sf (fiducial): yinit(H2) = 5 x fO^^ (solid), S5: 
2;i„it(H2) = 2 X fO-3 (long-dashed), S6: yinit(H2) = IQ-^ (short-dashed) and S7: yinit(H2) = IQ-" (dotted) 
are shown. 
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Fig. 6. Luminosity evolution of strong H2 and HD lines, (a) Those of the two H2 liens, 1-0 Q(l) 
(thick lines) and 0-0S(3) (thin lines), are shown for models SI (fiducial; solid), S6 (dashed) and S7 
(dotted), (b) Those of the same H2 lines (solid and long-dashed lines) as (a), and of the two HD lines, 0-0 
R(l) (short-dashed) and 0-0 R(7) (dotted) are shown for model S6. Formation of the central protostar 
corresponds to i = 0. Evolution up to the protostar mass of 70A/q is shown. 
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Evolution of the filamentary cloud of model Fl, where the contraction is induced by the H2 
Shown are (a) the temperature, (b) the H (solid) and D (dashed) molecular fractions and (c) the 
rate by H2 line (solid) and HD line (dashed) at the center. 
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Fig. 8. Evolution of the filamentary cloud of model F2, where the contraction is induced by the HD 
cooling. Meaning of the panels is the same as figure 7. 
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Fig. 9. Evolution of the filamentary cloud of model F3, the case of a high-density filament. Meaning of 
the panels is the same as figure 7. 
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Fig. 10. Luminosities of individual H2 lines (solid lines) and HD lines (dotted lines) averaged over the 
contraction time. Shown are for models (a) Fl, (b) F2, and (c) F3. The mass is normalized to lO^M© 
for all cases. Also, in panel (d), the case of a sph(2'S:al core in the collapse phase with the same initial 
condition as model F3 is shown. 
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Fig. 11. The luminosity in the H2 lines (solid) for a spherical core (the fiducial model SI) in the collapse 
phase against the number density. Dotted and thin lines show the values by the analytical expression 
(equation 29) for the temperature evolution in model SI and for constant temperatures of T = 500K and 
lOOOK. 
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Fig. 12. Total luminosity for filaments of model Fl (long-dashed; a low-density filament with efficient H2 
cooling), F2 (short-dashed; a low-density filament with efficient HD cooling), F3 (dotted; a high-density 
filament). Also shown by solid lines arc the analytical values (equation 36) for T = 100K,500K, and lOOOK. 
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Fig. 13. The redshift-star formation rate (SFR) domain wliere the strongest H2 hne flux F = Lci/47rL'J_|_^ 
exceeds the detection limit of SPICA, ~ lO^^^Wm^^. Tlie dark (hght) gray regions are the do- 
main detectable by SPICA for the cases where the luminosity from the whole pre-galactic cloud is 
Lci = lO'^^i/'iiierg s^^ (= 5 x lO'^'^V'iiisrg s~^, respectively). Pop III star-forming galaxies at redshift 
^ = 3 - 10 with Pop III SFR ~ lO^-^Mgyr"! are detectable by SPICA. 
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